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Abstract. Dual scattering channel (DSC) schemes generalize Johns' 
TLM algorithm in replacing transmission lines with abstract scattering 
channels in terms of paired distributions. A well known merit of TLM 
schemes is unconditional stability, a property that is commonly drawn 
upon the passivity of linear transmission line networks. So the question 
arises, if DSC algorithms remain stable in a neat sense. It is shown 
that a large class of a-passive processes are in fact unconditionally sta- 
ble. The analysis applies to TLM and DSC schemes alike and includes 
non-linear situations. 
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1. Introduction 

Dual scattering channel (DSC) schemes result from an incisive revision of 
the transmission line matrix (TLM) numerical method. The latter has origi- 
nally been introduced by P.B.Johns and coworkers in the early 1970s |JoBej 
and has since been subject to assiduous study and publication [Tlml-3]. The 
TLM method is today commonplace in scientific computing and largely ap- 
plied to the numerical solution of Maxwell's equations |Hoej but also to 
manifold wave propagation, transport, and diffusion phenomena; we can 
here refer to the monographs of Christopoulos jChj and de Cogan |dC| . 
Also, Rebel jRej gives a fairly complete survey over the state of the art of 
TLM by the year 2000. 

DSC schemes are generalized TLM methods in arising from a twofold 
abstraction |Hel| , Firstly, the scattering channel concept underlying TLM 
is redefined in terms of paired distributions. (Characteristic impedances 
are thus neither needed, nor in general defined, e.g.) In the second place, 
non trivial cell interface scattering is admitted during the connection step of 
iteration, thus taking advantage of the intrinsic duality in the connection- 
reflection cycle of the algorithm. 

The extended framework bypasses a set of modeling limitations induced 
by transmission lines ( and discussed in more detail in |Hel| section2] ) while 
it preserves the main advantages of the TLM method. In particular, the 
convolution type updating scheme and Johns' two step connection-reflection 
cycle are essentially retained. The question arises, however, if DSC schemes 
remain unconditionally stable in a well-defined sense. 

TLM schemes are unconditionally stable in that they are equivalent to 
passive linear transmission line network models |Jo2j . A concept that applies 
to DSC and TLM algorithms alike is a-passivity which paraphrases con- 
traction properties with respect to a non-negative (de)limiting functional; 
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cf. section |21 In the TLM context, a is essentially the sum over the squared 
incident and outgoing transmission line voltages. The sum is contracted at 
each scattering event due to energy conservation (or loss). 

In a like manner a-passivity characterizes the reflection and connection 
maps of an unconditionally stable DSC algorithm, but the 'energy' func- 
tional a (which may in fact measure any conserved quantity) needs not to 
be a quadratic form. For instance, energy may be linearly related to tem- 
perature and quadratically to particle velocity, within the same algorithm. 
The following approach is sufficiently general to apply to such and other, 
even non-linear situations. 

2. Stability 

Algorithm stability prevents the computational process from piling up to 
infinity (it does not yet imply convergence, or consistence of the algorithm, of 
course). TLM models are unconditionally stable in that they are equivalent 
to passive linear transmission line networks |.Io2j . DSC schemes, in not using 
lines, need a more general characterization which here is given in terms of 
a-passive causal functions. 

Let & be a real or complex linear space and / a totally ordered set ( e.g. 
/ € {N , Z , R} ; intervals are then naturally defined in I by the order re- 
lation. We commonly think with / of a discrete or continuous time do- 
main.) Also, let £ C £/ be a set of functions such that / £ £ implies 
X s<t ( s )f( s ) ^ f° r every t £ I , where X P ( S ) denotes the characteristic 
function of property P (which is 1 if s shares that property and else). 

Definition 2.1. A function F : £ — ► £ is called causal, iff for every / £ £ 
and t € I 

Ff(t) = F[ Xs<t f(s)](t) . 
Such functions are also called (causal) propagators. 

Remark. In some respect, causal functions generalize lower triangular ma- 
trices or integral operators such as 

Ff(t) = I K(t - s)f(s)ds 

J — oo 

with a Green's function kernel K, e.g. Note that in general F : / i— > Ff 
needs not to be linear. 

Typically / represents a state evolving in time (i.e. a process). Then 
causality of F means that F f (t) depends on the history of / only up to 
present time t . 

The proof of the following is easy and left to the reader. 
Proposition. 

For every t € / ; /, g € £ and causal functions F, G : £ — > £ 

(i) F f ( s ) = F [ X rsit (r) f (r)](s) , for every s < t . 

(ii) V f(s) = g(s) for every s ^ t , 

then Ff (s) = Fg (s) for every s ^ t . 
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(iii) The product of causal functions 
FG : £ — > £ 

/ FG f := F[G[f}) 

is again causal. In fact, if £ is a linear space, then the causal func- 
tions over £ form an algebra. 

Let ||...|| be a norm on XL and a G M.^ a continuous non- negative real 
functional on XL. 

Definition 2.2. 

(i) We call a process g : 7 — ► XL stable, iff 5 is uniformly bounded on 7 , 
i.e. iff there exists b G M+ such that || g (t) \\ < b for every t G 7 . 

(ii) The functional a : XL — > M is named a (de)limiting functional, iff 
there exist any non-negative real constants a, b, c such that 

(1) ||z|| < a + 6(a(z)) c , 

for every z G XL . Then obviously 6, c > , and we say also that a 
is minimal increasing {in any order) not lesser than 1/c. 

Let pi be a measure on 7 such that intervals are /i-measurable sets. 
Functions on 7 are henceforth read modulo fi (viz. as equivalence classes of 
functions that differ at most on sets of /x- measure zero). Also, let ael 1 
be a delimiting functional on XL (i.e. one that is increasing not lesser than 
any positive order 1/c), and assume that a o f is jjl— summable over finite 
intervals in 7 for every / G £ . The latter is for instance the case if a(z) = 
|| z || p for any real p ^ 1 and £ C L p ( 7, XL ) , which is the metric completion 

° f { / e XL 7 I (J\\f\\ p d»)* < 00} , 

i.e. £ is a subset of the Banach space with norm || / || p := ( f T || / || p dfi ) p . 
Definition 2.3. A causal function 7 1 : £ — > £ is called a-passive, iff 

(2) / a(F/(s))d/x(a) < / a (/ (« ) ) d^W > 
for every / G £ and f 6 7. 

Remark. If a = ||...|| p for any real p ^ 1 and || / || p is /i-summable over 
7 , i.e. / G L P (I, XL ) , then @ clearly implies 7 1 / G 7 P ( 7, XL ) and 

11*711, < II /lip • 
Hence, every ||...|| p -passive causal function F defines a contraction operator 
on £ n 7 p (7, XL). 

Assume, furthermore, that r G M+ and let on 7 : = M the measure (i be 
concentrated in {kr | A; G Z } with uniform weight /i ( {/cr} ) = r , k G Z . 
Alternatively, let 7 : = { /cr | k G Z } with 'the same' measure . 
(Virtually we deal of course with that discrete situation, even in working on 
the real axis with functions that are constant over intervals [kr, (k + 1)t) ; 
k G Z . Indeed we retain the integral formalism for simplicity, and the reader 
may optionally re-write the following integrals as sums.) For every / G £ let 
f(t — t) G £, i.e. £ is closed under time shifts by negative integer multiples 
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of r. Then, for arbitrary N E N+ and exciting function e E £ C £ with 
support on [0, A^r ) C /, the following holds 

Theorem 2.1. [Stability of the iterated passive causal process] 
For any a-passive causal function F : £ — > £ and e E £ as stated, 
if g E £/ is a process such that g(t) = /or i ^ and recursively for 
t = nr; n E N 

(3) 5 (t +r) = F[e + <?](i) , 

i/ien (7 is uniquely defined (modulo fi) and for every t ^ Nr holds 

b f 

(4) < a + ( - / a(e + g)-a(g) dfi ) c 

T J[Q,Nt) 

with every constants a,b,c that satisfy ©. Hence g is stable. 

Remark. A process g E £/ which is recursively generated according to © 
by iteration of an a-passive causal function F is called an a-passive process. 

The theorem ensures thus that for any excitation of finite duration (and 
with no further restrictions) the a-passive process is necessarily stable and 
in this sense unconditionally stable. 

Note that existence of such a process g is not apriori guaranteed, since 
this obviously depends on the condition that with so : = e E £ recursively 
also the functions 

(5) s n (t) : = e(t) + F[s n -i ](t — t) E £ , 

for < n < N - which eventually has to be checked. 
Clearly, conditions © are always true (and hence g exists) if £ is a linear 
space. We do not universally premise this, in order to potentially apply the 
theorem to non-linear situations, where conditions © may only be satisfied 
for sufficiently small excitations e . 

Corollary. 

(i) In the special case a = ||...|| estimates Q holds with a = 0, b = 1, 
c = 1 . Then the triangle inequality applies to the integrand of (@J) 
and validates the bound 

\\g(t) || ^ - / ||e|| dpb . 

T J[0,Nt) 

(ii) If N = 1 , i.e. e(t) is a DlRAC excitation concentrated on [0,r) 
(where g = 0), then @ reads simply 

\\g(t)\\ ^ a + ( - / a(e) d\i ) c , 

provided that a ( ) = ( which is the normal case ) . 

Proof. Clearly, g E £ is uniquely defined by the given recurrence relations, 
since e + g ( and hence g ) at the right hand side of © is evaluated only 
up to time t = n r , in virtue of the causality of F . 

Furthermore, if N ^ n and a ^ \\g (nr) || with any a satisfying JQ), then 
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with pertinent b, c that satisfy (JJ) 

< (U\\g(nr)\\ - a)) 1 ^ ^ a(g(nr)) 



b 



= -( J a(g)dfi(s) - I a(g)dfi{s)) 

T Js<(n+1)T J s<nr 

= - ( I a(F[e + g}) dfi(s) - f a(g) dfi(s) 

1~ J s<nr J s<nr 

< recursion formula @ 

^ - ( / a(e + g) dfj,(s) - j a(g)dfi(s)) 

s<nr J s<nr 



a(e + g) - a{g) d(i(s) 

0,Nt) 



T 

t since F is passive 
1 

T 

t since e(t) = if t <£ [0, Nt) . 

Thus, estimates (@J) holds true in the case a ^ ||</(nr)|| and trivially 
otherwise. It follows that || g || is uniformly bounded on / \ [0, Nt) , hence 
also on / , that is to say g is stable. □ 

3. DSC Processes 

In this section DSC schemes are represented as paired a-passive processes 
such as are dealt with in Theorem 12.11 and which hence are stable. 

Just as the TLM algorithm, DSC schemes operate on a space T of prop- 
agating fields, which is a product of (real or complex) normed linear spaces 

(6) T = Tin X T Q ut i 

cf. |Helj . The two factors are named (somewhat off-hand) the incident and 
outgoing subspace of T. They are isomorphic in that there is canonical 
involutary isomorphism of normed linear space 

nb : T -» T 

(7) 

Z — ( %in j Zout ) 1 * ( %out 5 Zin ) — • ^b ( Z ) , 

which is commonly called the node-boundary map. Hence, there exists a 
space (£j, ||...||) such that 

(8) Tin — T ou t — II---II) 

in the sense of isomorpy of normed spaces and 

? = (£ 2 , ll-ID , 



e.g. with norm || ( a , b) ||~ : = \/ll a ll 2 + ||6|| 2 ; a,b £ L 

(or any equivalent norm). 

As is well known, DSC and TLM algorithms follow a two-step iteration 
cycle in working with alternate application of a connection and reflection 
map 

G : I X T ou j. > P-in and 3£ : ,7 x IPj n ► Tom* j 
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which respectively update the propagating fields at even and odd inte- 
ger multiples of half the time step, i.e. on J := {kr \ k £ Z } and J : = 
{ (2k + 1)t/2 I k G Z }. To these maps the following functions F c , F R are 
associated in a one-to-one correspondence 



(10) 
and 
(11) 



F„ : L 1 -> L 1 



c 



f » F c f with F c f(t) := nboe(t, f ] 



F R : L J -» £ 



g ^ F R g with F R g(t) : = nbo%(t, g) 



Proposition 3.1. 

(i) -F c and i 71 ^ are causal on £/ and £ , respectively. 

(ii) For r, s £ J and T s : /(£) i— ► f(t + s) the shift operator on £ 



I u J 




F 

( r, s ^ and < c is a-passive on 

, \ T r o F„ oT s . 
=> ( < „ „ 2S a-passive on 

Indeed, the first statement is only a trivial consequence of the definition 
of the DSC propagators 3? and C as functions on back-in-time running 
sequences of incident and outgoing states, cf. |Helj . section 3. Then (ii) 
follows from Definition 12.31 □ 

Definition 3.1. 

„, connection map C . „ , .„. I F n 

Ine < n . m is called a-passive, iff < ^ is a-passive 

I reflection map Ji b 



in the sense of Definition 127 
We now claim the main statement: 

Theorem 3.1. With every time limited excitation, the DSC process gen- 
erated by a-passive reflection and connection maps is uniformly bounded, 
hence stable. 

Proof. It is sufficient to show that every finitely excited DSC process that 
is generated by a-passive 31 and C can be written as a pair of processes, 
either of which satisfies Theorem 12. II 

With H : = I U J = {kr/2 \ k € Z} and the measure on H inherited 
from / and J jointly, the space of all DSC processes is 

£ := {h=(f,g)e(L 2 ) H \ f((2k + l)r/2) = f(kr) 

and g(kr) = g((2k — 1) r/2 ) , for every k £ Z } , 

i.e. the functions / : H — ► £ and g : H — > £ in h = (f,g) £ £ 'switch' 
at even and odd integer multiples of r/2 , respectively. So, there is a natural 
bijection £ ^ £j x lJ 

h = (f,g) ^ (f{I,g[J) , 
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in virtue of which the first and second components in £ can be naturally 
identified with L and & J , respectively. 

For any incident function e G £ x { } C £ supported on a finite 
interval [ , Nt ) C I , the DSC process excited with e and generated by 
C and 01 is the well-defined function h G (L 2 ) H such that fo(t) = for 
i ^ and recursively for < t G if 

(12) fc(t + l ) _/<M*).r- i *.T_ i [ e + * 1 ](t)) if«s/ 

Actually h is well and uniquely defined by relations 1)12(1 which provide 
separate recurrence relations for the two processes hi , h 2 ; for instance for 
h\ and t G I 

h 1 (t + t ) = /ii ( r + - ) with t~ : = t + ^ G J 



T=rF c T=r[h 2 \ r ](i~) by© 

T , 



2 C 2 



(13) = ^[J^Kt 

= TrMj- . cf. G2J 
= F ZU.F H :ZU[/n + e] (t) . 

In virtue of the causality of F C T- T F R T-t function hi enters the last 

expression in (|13|) only up to argument t . Moreover, since products of 
a-passive operators are obviously a-passive, F c T_ T / 2 F R T_ T / 2 is a-passive 
due to Proposition 13. II (ii) . Hence Theorem 12.11 applies to h\ , which thus is 
stable, just as then also is h 2 = F R T_ T / 2 [e + hi]. □ 

4. Conclusions 

It has been demonstrated that a class of a-passive processes, charac- 
terised by simple contraction properties with respect to a limiting functional 
a , are unconditionally stable. Finitely excited DSC processes, generated by 
a-passive reflection and connection maps, are then necessarily stable. 
Stability analysis reduces so in pursuant cases to finding any limiting func- 
tional a in relation to which the reflection and connection maps are 
a-passive. In reverse: If the connection and reflection maps are per se 
so designed as to conserve ( or admit loss of ) any quantity 'measured' by 
a limiting functional, then the DSC process generated by those maps is 
essentially stable. 
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